Electrostatically induced phase transitions in superconducting complex oxides 
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We describe quantum phase transitions in superconducting complex oxides which could be tuned 
by electrostatic charge transler. Using a simple model for the superconductivity of a thin film or 
surface of a bulk copper oxide, we show that tuning the carrier density may allow the visitation 
of several superconducting phases with different pairing symmetries such as extended s- (se), d- 
and (se ± id)-wave. We construct a universal phase diagram for single-band superconductors with 
se- and d-wave components of the order parameter based on symmetry considerations alone. For 
a specific model with nearest neighbor attraction, we obtain the phase diagram in the interaction 
versus filling factor space showing the boundaries of the possible phases. Finally, we calculate the 
superfluid density and penetration depth as characteristic properties of each phase. 
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PACS numbers: 74.78.-w, 74.78.Bz, 74.62.-c, 74.62.Dh 

Substantial progress in tuning carrier density has been 
achieved recently in systems like amorphous Bismuth, 
where a superconducting-insulator transition was driven 
by electrostatic fields 1 . This experimental realization has 
opened the possibility of studying quantum phase tran- 
sitions by tuning carrier density electrostatically. Sev- 
eral groups are currently attempting the use of electro- 
static doping to control the carrier density in complex 
oxides, with particular attention to cuprate superconduc- 
tors. The phase diagram of cuprate superconductors as 
a function of chemical substitution (oxygen deficiency) is 
known for Lanthanum (Yttrium) cuprates, but as it is 
well documented in the literature that the relation be- 
tween chemical doping or oxygen deficiency to carrier 
density is very complicated 2 . 

Eletrostatic doping has substantial advantages over 
chemical doping, as the carrier density can be tuned con- 
tinuously and does not produce undesired effects such 
as disorder or changes in the crystal structure as chemi- 
cal doping or oxygen deficiency often does. The method 
when applied to cuprates should be able to answer the 
long standing question of what really is the phase dia- 
gram of the system as a function of carrier density, and 
to clarify the origin of the dome that reflects an optimally 
doped superconductor with the highest critical tempera- 
ture. In particular, if the results of electrostatic doping 
are substantially different from those of chemical dop- 
ing, showing for instance that the critical temperature 
continues to grow beyond the expected optimal chemi- 
cal doping point, then it could be directly inferred that 
chemical doping does introduce undesired effects. 

Although electrostatic doping has remarkable advan- 
tages over chemical doping, its use is limited to thin 
films and the surface of bulk systems. However, it is 
probably the best way to study the dependence of phys- 
ical properties of complex oxides and other materials on 
carrier density. Furthermore, from the point of view of 
devices, the electrostatic tuning of the phase diagram at 
fixed temperature would be extremely useful to produce 
a superconducting switch, as the resistance of a thin film 
could be turned on and off via a gate voltage. 



It is in anticipation of these exciting experiments cur- 
rently underway in several labs around the world, that 
we describe the possibility of using electrostatic tuning 
of carrier density to study quantum phase transitions in 
complex oxides. As an example we chose to study a sim- 
ple model of complex oxides representing a cuprate su- 
perconductor. We model such a system via an extended 
attractive Hubbard model in a two-dimensional lattice, 
and derive the universal phase diagram when only two 
pure order parameter symmetries are allowed. Based on 
symmetry alone there are several possibilities: pure ex- 
tended s- (se) and d-wave phases, and (se ± d) phases 
that do not break time-reversal symmetry and (se ± id) 
phases that do. For a specific model where only nearest 
neighbor attraction is included, the se-wave phase dom- 
inates at lower filling factor, while the d-wave solution 
dominates at higher filling factor, with a (se ± id)-wetve 
phase in between, whereas the (se ± d)-wave is not ac- 
cessible. The phase diagram for such a minimal model is 
very rich, and the existence of various phases as a func- 
tion of filling factor for fixed interactions can be tested 
via measurements of the penetration depth (superfluid 
density) at low temperatures. 

To study the physics discussed above we start from a 
two-dimensional hamiltonian 
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describing thin films or the surface layer of a bulk cuprate 
on a square lattice. Here, the local particle number op- 
erator is rii a — ctj.Cjo-, while — U and — V correspond to 
on-site and nearest neighbor attractions, respectively. 

In order to establish the quantum phases as a function 
of carrier density, we start by constructing the partition 
function Z = J "Dc'Vce for the action 
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dr 
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Introducing order parameters for local s-wave A s / , ex- 
tended s-wave A se and d-wave A^ pairing, we can write 
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the action as 
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where a = sl,se,d; the interactions are T4; = U, V st 
V c i = V; and the four- vector q — (ii/ n ,q). 

The inverse free fermion propagator matrix is 
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with kinetic energy £k = £k — /x, band dispersion 
6k = — 2i [cos(fc x a) + cos{k y a)], chemical potential \x, 
four- vector k = (zw„,k), and unit cell length a. The 
additional matrix appearing in Eq. ([3]) is 
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where the Einstein summation over a is understood, and 
A Q (k, k') are the symmetry factors for the order param- 
eters, which in the limit of zero momentum pairing (k = 
k') become A s ;(k, k) = 1, A se (k, k) = cos(fca;a)-|-cos(fc a a), 
Ad(k, k) = cos(fc x a) — cos(k y a). 

The saddle point approximation gives qualitatively 
correct results either at low temperatures (T <C T c ) 
for any interaction strength or for weak interactions 
V <C bw = 8t at any temperature^, where bw is the 
bandwidth in two-dimensions. In these cases, the effec- 
tive action can be approximated by 

S = -^|A si | 2 - ^(|A se | 2 + |A d | 2 ) + S 2 + ^. (6) 

The second term in the action 52 = 
Sk 7 m [1 + ex P(~ Ektf/T)], contains the quasi- 
particle (7 = 2) or quasihole (7 = 1) energies 



i?k, 7 = (— ) 7 y ^k + |A Q A Q (k)| 2 , where the symmetry 
function is defined as A Q (k) = A Q (k, k). 

Notice that even in such a simple model with local 
(-U) and nearest neighbor (— V) attractions the num- 
ber of possible phases is quite large. There are three 
possible pure phases: local s-wave {si); extended s-wave 
(se) and d-wave (d). There are several possible binary 
mixed phases si ± se, si ± d, and se ± d, which do not 
break time-reversal symmetry, as well as those that do 
like si ± ise, si ± id, and se ± id. Lastly there are also 
several ternary mixed phases involving all three symme- 
tries si, se, and d. Since the situation is quite compli- 
cated in the more general case, for clarity and simplicity, 
we study here the case of U = 0, which allows for solu- 
tions involving only se- and <i-wave symmetries. In this 
simpler case there are only six possible phases. The pure 
phases are se and d. The mixed phases are se ± d, and 
se ± id. Thus, from now on, we will confine ourselves to 
this simpler analysis. In this case, the order parameter 
equations can be obtained by minimizing the action S 
through the conditions SS/SA* se = and SS/SA* d = 0. 



For the se-wave component of the order parameter the 
first condition leads to 
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where A se (k) = A se [A se (k)] 2 + A d A se (k)A rf (k). Corre- 
spondingly for the d-wave component 
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where A d (k) = A d [X d (k)} 2 + A se X se (k)X d (k). 

The number equation that fixes the chemical poten- 
tial is obtained through the thermodynamic relation 
N = —dft/dfi, where ft = —T\nZ is the thermodynamic 
potential. In the present approximation the thermody- 
namic potential is Q = —TS, and the number equation 
reduces to 
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where v = N/N s is the filling factor. 

The phase difference between the s-wave and d-wave 
components of the order parameter is defined to be 
S<f> = 4>d — 4>se, where <pd is the phase of the ei-wave or- 
der parameter A d = \Ad\e l ^ d and <j> se is the phase of the 
s-wave order parameter A se = |A se |e*^ se . Simultaneous 
solutions of Eqs. (J7J), ©, and © reduce to saddle-point 
solutions of only for 5(f) = tt/2, 3tt/2 which correspond 
to se ± id phases that break time-reversal symmetry, and 
Sep = 0, 7r which correspond to se ± d phases that do not. 

The saddle point critical temperature can be obtained 
by setting the order parameters A se = and A^ = 
in Eqs. (O, {SJ, and ([9]). In this case, the filling fac- 
tor dependence of critical temperature T c (y) and critical 
chemical potential nd^) can be obtained for pure se- and 
d-wave symmetries. The solutions for T c {y) are shown in 
Fig. Q] for V/t = 3.0. Notice that the s-wave phase is 
favored at lower filling factors, while the e?-wave phase 
is favored at higher filling factors. In addition, the crit- 
ical temperature is symmetric about v — 1, since the 
Helmholtz free energy F = £1 + fJ,N is invariant under 
the transformation fj, —/j, and v —> 2 — v. 

In order to construct a universal phase diagram for all 
possible phases it is important to construct the Ginzburg- 
Landau theory near T c by expanding the action described 
in Eq. © in terms of the order parameters A se , A^ and 
their complex conjugates. From the thermodynamic po- 
tential fl = —TS, we can calculate the Helmholtz free 
energy F = O + fj,N. The free energy per site T — F/N s 
takes the simple form 



J~ ^se|A se 



- a d \A d \ 
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This expression has precisely the expected form, based 
on symmetry grounds alone, when paring with se- and 
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FIG. 1: Plots of the critical temperature T c as a function of 
filling factor v at fixed interaction (V/t = 3.0), showing the 
the se-, d-, and (se-Md)-wave phases. Notice the tetracritical 
point where the normal and all superconducting phases meet. 



d-wave components are possible. However, the coeffi- 
cients a and b depend explicitly on the parameters of the 
microscopic model used. For the construction of the uni- 
versal phase diagram we introduce the dimensionless pa- 
rameters X — {b s d [l + 5 cos(2Scf))\ /bd} x \a,d/a se \, and 
Y = (b se /bd) x \ad/a se \ 2 . Minimization of the free energy 
in Eq. (|10j) and a stability analysis leads to the universal 
phase diagram shown in Fig. [2j where all possible phases 
(se, d, se ± d and se ± id) are indicated. Notice that the 
free energy depends only on 28<j) and does not distinguish 
between the phases se + d and se — d, which are thus de- 
generate. The same applies to the phases se + id and 
se — id, which are also degenerate. In the particular case 
of se±id phases, time-reversal symmetry is broken while 
chirality is not, the latter of which requires additional 
terms in the free energy for the distinction between the 
se + id and se — id phases. Notice that a tetracritical 
point occurs at Y = X = 1, where the normal and su- 
perconducting phases with se, d and se ± id symmetries 
meet. 
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FIG. 2: Universal phase diagram for superconductors ex- 
hibiting order parameters with se- and d-wave components is 
shown, in the space of the dimensionless parameters X = 
{bsd [l + | cos(2<5</>)] /b d } x \a d /a ae \ and Y = (b se /b d ) x 
\ad/a S e\ 2 - Depending on microscopic details, different regions 
of the phase diagram are accessible. 

For the specific microscopic model discussed above the 



variable X is always positive, such that the system is 
always stable and the se ± d phases are not accessible. 
Thus, we elaborate further only on the accessible phases 
se, d, and se ± id. At low temperatures T <C T c , the 
Helmholtz free energy per site can be approximated by 

? = \ I A «! 2 -^ E ln I 1 + exp(-S k , 7 /T)]+ M (^-l) 

a k,7 

for any interaction strength V/t. A comparison of the 
free energies for all the accessible phases produces the 
phase diagram in the interaction (V/t) versus filling fac- 
tor (v) space shown in Fig. [3] Notice that pure se and 
d phases are always separated by (se ± id) phases, and 
their phase boundaries describe continuous transitions. 
For weak attractions the region of filling factors where 
(se ± id) phase is realized is very narrow. However, the 
region increases substantially as the interaction V/t gets 
larger. For instance when V/t = 3.0, the (se ± id) phase 
exists between v m i n w 0.36 and v max w 0.48. 

We had hoped that a topological quantum phase tran- 
sition^ characterized by the emergence of a gapfull d- 
wave superconductor from a gapless d-wave superconduc- 
tor would also emerge within the d-wave region of the 
phase diagram. However, within the d-wave boundary 
the chemical potential always fall inside the band limits 
(|/i| < 4i), and one can always find zeros of the quasipar- 
ticle excitation spectrum given by £k,2 = 0. However, 
the transition from d- to (se±id)-wave is also very exotic 
as it involves a change in the excitation spectrum from 
gapless to fully gapped with a corresponding change in 
topology of the quasiparticle-quasihole excitation mani- 
fold, and a change in order parameter symmetry accom- 
panied by the spontaneous breaking of time-reversal. 

Since in standard condensed matter physics the inter- 
actions in the same material can not be tunned, one can 
hope to visit different phases by changing the tempera- 
ture or tuning the filling factor (carrier density)^. If the 
tuning of carrier density via electrostatic means can be 
achieved experimentally for complex oxides, then quan- 
tum phase transitions may be studied as a function of fill- 
ing factor—. Since electrostatic tuning of carrier density 
can only be implemented in thin films or at the surface 
of bulk materials, additional experiments to characterize 
the various phases are difficult. However, in such ge- 
ometry, measurements of the penetration depth \(y, T) 
may be possible. Given that \{y,T) cx p~ 2 {v : T), where 
p{v, T) is the superfluid density, then a distinguished low- 
temperature behavior varying from a linear increase with 
temperature in the d-wave phase to an exponentially ac- 
tivated behavior in the se-wave phase could be revealed 
as the filling factor v is varied. 

In order to describe the behavior of \{y, T), we calcu- 
late the superfluid density tensor 

Pij (y, r) = ^E [2n k 0i0i& - H'UVUk , (11) 

k 

at low T, where n k = (1/2) [1 - ^ta,nh(E k , 2 /2T) / E K2 ] 
is the momentum distribution, and Y"k = 
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FIG. 3: The zero-temperature phase diagram in the interac- 
tion V/t versus filling factor v plane showing the boundaries 
between all accessible phases se, d, and se + id for the mi- 
croscopic model discussed. The phase transitions are always 
continuous across the boundaries. 

(2T)" 1 scch 2 ( J B ki2 /2r) is the Yoshida function. Notice 
that p xx = p yy and p xy = p yx for all order parameter 
symmetries. In Fig. [4j we show the temperature depen- 
dence of the superfluid density in the se-, (se ± id)-, 
and d-wave phases. At low temperatures, the superfluid 
density for the se-wave phase exhibits an exponentially 
activated behavior p se {0) — p se (T) ~ exp(— |A se |/T), 
due to the presence of full gap in the quasiparticle 
excitation spectrum i?k,2, with similar behavior for 
the (se ± id) phase. However, in the d-wave case the 
superfluid density decreases linearly with temperature 
Pd(Q) — Pd{T) ~ T, as expected from the nodal structure 
of£ k , 2 . 
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FIG. 4: The superfluid density p/t versus temperature T/t 
at fixed interaction V/t = 3.0. The dotted line corresponds 
to se-wave for v = 0.24, the solid line to (se ± id)-wave for 
v — 0.38, and the dashed line to d-wave for v — 0.60. 



Before concluding, we would like to make an important 
remark. For the model described here there is no antifer- 
romagnetic phase near half-filling. This phase emerges by 
replacing the on-site attractive-?/ term by an on-site re- 
pulsive (Hubbard-t/) term in the Hamiltonian of Eq. fl} . 
The full solution of the problem including the on-site 
repulsion U, and the nearest neighbor attraction (— V) 
is quite complex, however one can make a few qualita- 
tive statements for U/t ^> V/t and fixed V/t. In this 
case, the system is an antiferromagnetic insulator at and 
near half-filling (y = 1), however away from it the ef- 
fects of a locally repulsive term are dramatically reduced 
and with decreasing filling factor superconductivity is 
achieved first for the d-wave phase, then for the mixed 
phase (se ± id) and finally for the se-wave phase. Since 
the change in symmetry of the order parameter from d- to 
mixed phase to se-wave occurs reasonably far away from 
half-filling, such a transition is not dramatically affected 
by U . Thus, even in more realistic models for complex 
oxides such as the cuprates, the transition proposed here 
should persist at lower filling factors. 

We have discussed a simple extended attractive Hub- 
bard model to describe single-band complex oxide su- 
perconductors in two-dimensions, where the filling factor 
can be adjusted via electrostatic doping. Based on sym- 
metry alone, we established that the possible supercon- 
ducting states are extended s-wave (se), d-wave (d) and 
mixed phases which break (se ± id) and do not break 
(se ± d) time-reversal symmetry. However, we found 
that only the se-, d- and (se ± id) phases are accessi- 
ble within a nearest neighbor attraction model, and that 
there exists a tetracritical point where the normal and 
all superconducting phases meet at finite temperature. 
We have shown that quantum phase transitions between 
various superconducting phases take place at filling fac- 
tors far from half-filling, and analysed the temperature 
dependence of the superfluid density (penetration depth) 
near the boundaries of such transitions, where the char- 
acteristic power law behavior of the d-wave symmetry is 
replaced by the exponentially activated behavior of the 
se-wave symmetry. 
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